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ABSTRACT 

Quasi-toroidal oscillations in slowly rotating stars are examined in the framework of 
general relativity. The oscillation frequency to first order of the rotation rate is not a 
single value even for uniform rotation unlike the Newtonian case. All the oscillation 
frequencies of the r-modes are purely neutral and form a continuous spectrum limited 
to a certain range. The allowed frequencies are determined by the resonance condition 
between the perturbation and background mean flow. The resonant frequency varies 
with the radius according to general relativistic dragging effect. 
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1 INTRODUCTION 

In recent X-ray observation with the Rossi X-Ray Timing 
Explorer (RXTE), quasi-periodic oscillations are discovered 
in several sources (e.g., van der Klis et al. 1996 and the sub- 
sequent papers of the volume). The frequency ranges from 
a few Hz to kHz and may be attributed to the phenomena 
near a compact object. Several models are proposed as the 
oscillations: beat-frequency between a magnetized neutron 
star and accretion disc, stellar oscillation, and so on. For ex- 
ample, Strohmayer and Lee (1996) considered the excitation 
of the g- and r-modes as a result of the thermonuclear flash 
and discussed the observational possibility. Their calcula- 
tions are however based on the Newtonian gravity. Unlike 
the spheroidal modes like f-, p-, and g-modes, the general 
relativistic effects are not clear for the r-mode, since the 
problem has never been studied so far. The toroidal motion 
is trivial in a non-rotating star, but has non-vanishing fre- 
quency in a rotating star. The quasi-toroidal mode is called 
r-mode and known as the Rossby wave in ocean. Papaloizou 
and Pringle ( 1978 ) introduced the r-mode in connection with 
the variable white dwarfs. See also the subsequent study 
(Provost, Berthomieu & Rocca 1981; Saio 1982). 

In this paper, we will explore the relativistic effect on 
the r-mode. We never discuss the observational implication 
of the r-mode in relativistic stars, but theoretical study of 
the oscillation frequency may be a useful tool for the future 
observation. We use the slow rotation approximation and 
linearized Einstein equations. The first-order effect of the 
rotation rate is taken into account. In section 2, we present 
the perturbation equations describing the r-mode. In section 
3, the eigenvalue problem is solved. Finally, section 4 is de- 



voted to the discussion. Throughout this paper we will use 
units of G = c = 1. 



2 PERTURBATION EQUATIONS 

We assume a star with a uniform angular velocity ^ 0(e), 
and consider the rotational effect of order e only. The con- 
figuration of the pressure p and the density p is the same as 
in the non-rotating star, since the centrifugal force deform- 
ing the shape is of the o rder e^. The metric for the slowly 
rotating star is given by (Hartle 1967) 



ds 



9^ + sin^ ( 



-2Lor'^ sin^ Odtdtj), 



(1) 



where u ^ 0{e) is a radial function describing the dragging 
of the inertial frame. Introducing a function vj = fl — ui, we 
have a differential equation as 

(^jr'^zu'^ ~ 16tv{p + p)e^jr'^-aj = 0, (2) 

where a prime means a derivative with respect to r, and 

j = e-^^+'-^Z^ (3) 

The function zu inside the star is monotonically increasing 
function of r, so that the range is limited to 



(4) 



where zuq and zur are the values at the center and surface 
{r = R), respectively. 

The perturbations describing non-radial oscillations 
with the small amplitude can be given by the density per- 
turbation Sp, pressure perturbation 5p, and three compo- 
nents of the velocity {U,R,V). The metric perturbations 
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can be expressed by the ten functions, but the number 
is reduced to six {ho, hi, Ho, Hi, H2, K) by the gauge fix- 
ing. We here use the same notation for these perturba- 
tion functions as in Kojima (1992), but the exphcit forms 
are not necessary for most of the foUowing discussion. In 
this way, the equations governing the osciUations are one 
thermo-dynamic relation and ten components of the lin- 
earized Einstein equations for these eleven functions. In the 
case of non-rotation, two sets are completely decoupled. One 
set (U,ho,h\) is called axial perturbation (or "odd-parity" 
mode), while the other set {5p, 5p, R, V, Ho, Hi, H2, K) polar 
perturbation (or "even-parity" mode). Notice that the axial 
perturbation U describes the toroidal motion, and has zero 
frequency in the non-rotating star (Thorne & Campolattro 
1967). We expect that with rotation the toroidal oscillations 
of the fluids have finite frequencies of the order of f2 like 
in the Newtonian pulsation theory (Papaloizou & Pringle 
1978; Provost, Berthomieu & Rocca 1981; Saio 1982). There 
also exists gravitational wave mode with non- vanishing fre- 
quency in the axial perturbation (Chandrasekhar & Ferrari 
1991; Kokkotas 1994). The mode can be distinguished from 
the r-mode in the non-rotating limit. We never discuss the 
gravitational wave mode any more here. 

We look for the r-mode oscillations in the relativistic 
rotating stars. The perturbation functions are expanded by 
appropriate sets of spherical harmonics with index I, m and 
exp[— i(crt — m0)]. The linearized Einstein equations in the 
slowly rotating star are schematically given by 



Aim + f A X Vl±lm ~ 0, 
Pim + X Al±lm = 0, 



(5) 
(6) 



where A and A represent some sets of the axial perturbation 
functions, while V a nd V represe nt those of the polar per- 
turbation functions (Kojima 1992). The symbols, £p and £a 
are some operators of the order e. It is clear that the presence 
of the rotation induces the couplings between the axial and 
polar modes. The coupling is subject to the selection rule: 
the axial mode with /,m is coupled with the polar modes 
with Z ± l,m and vice versa. This rule is easily understood 
if we notice that the slow rotation perturbation corresponds 
to the odd-parity perturbation with I = 1 (Campolattro & 
Thorne 1970). 

In the previous papers (Kojima 1992, 1993a,b), the pul- 
sation equations in the slowly rotating stars are examined, 
assuming that the oscillation frequencies in the non-rotating 
stars are regarded as non-zero values. This is true for the 
spheroidal modes like f-, p-, g- modes and gravitational wave 
modes. The eigenvalue problems are solved for the non- 
rotating stars, and the rotational corrections are calculated 
for these oscillation modes. We instead assume that the fre- 
quency a is of the order of e. Different manipulation is there- 
fore necessary, since the rotation should be included at the 
lowest order to obtain the oscillation frequency of the r- 
mode. The perturbation functions should be ordered in the 
magnitude as, 

-0{e^U), 

-0{e'U), 
'O(e-'U), 



ho ~ OiU), 
Spr^Oie^U), 
ffo ~0(e'(7), 
Hi ~ Oie^U), 



hi r 

H2 



K 

V ' 



'0{e^U), 
0{e^U). 



(7) 



The velocity perturbations and the metric perturbations 
with htj are 'anti-symmetric' with respect to time and oth- 



ers like density perturbations are 'symmetric'. The former 
should therefore have even power of e, while the latter odd 
power of e. The polar perturbation functions should be 
of higher order in the quasi-toroidal oscillations. From the 
above ordering, the linearized Einstein equations correspond 
to those for the axial part, Aim = at the lowest order. The 
polar part is induced through the coupling, equation (^) at 
higher order level. The corrections to the axial part are also 
induced at higher order level. In proceeding to the higher 
order, the higher order rotational corrections for the equi- 
librium states are necessary, but the corrections of the order 
e are sufficient at the lowest order. 

We now solve Aim ~ for U, ho, hi. The quasi-toroidal 
velocity can be expressed as 



a — mil + 



2mvj 



U = -4n{a ~mn){p + p)e"'ho. (8) 



The relation between the metric perturbations is 



hi = - 
where 



(9) 



(10) 



The master equation governing the quasi-toroidal oscilla- 
tions can be written as 



(ro - fi) 



where 



+ 1) - 2] , 

1 = -4 U'^ ^ ) 

= 16n{p + p)e^zu, 
and the eigenvalue 

-{a — mil). 



2m 



(11) 

(12) 

(13) 
(14) 

(15) 



In equation (0), we have used the relation 



3 SINGULAR EIGENVALUE PROBLEM 

The basic equation ( |ll| ) is not a regular eigenvalue problem. 
The coefficient {vj — /i) becomes singular inside the star for 
a certain value of p.. The coefficient also vanishes outside the 
star, but the singularity can be removable because g = 0. 
This equation is very analogous to the Rayleigh's equation 
for the incompressible shear fiow (e.g., Lin 1955). The per- 
turbation propagating with the wave number k and speed c 
in the mean flow with velocity it can be described as 



(16) 



The similar singular eigenvalue problems appear in many 
other flelds, e.g., differential rotating ffuid dis cs and plasma 
oscillations. See e.g., Balmforth & Morrison ( 1995 ) for the 
methods of solving the singular eigenvalue problem. The sin- 
gular point is called as critical layer in the fluid dynamics, 
or co-rotating point in the rotating discs. The studies of 
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the singular eigenvalue problem indicate that unless there 
is an inflection point, u" = 0, somewhere within the flow, 
the eigenvalue is not discrete, but continuous and neutral 
against the stability. All neutral modes must have critical 
layers (co-rotation points) that lie within the flow, and there- 
fore form a continuous spectrum of intrinsically irregular 
eigen- functions. 

The parallel argument holds for our problem. The es- 
sential points are that the potential v is positive deflnite for 
I > 2 and that there is no inflection point {q > 0) inside the 
stars. We can conclude that the eigenvalue of equation jlll ) 
is real number and the range is limited to 



Tuo < jj, = (a — mil) < vjR, 

2m 

where the range of zu, M) is used. 



(17) 



We shall simply show the conclusion by reductio ad ab- 
surdum. If there is a non-trivial solution of which the eigen- 
value ^ is not located within the domain ( |l7| ) , then we have 
the integral relation 







-'- l/F.|2 • 4 , 

jr ar, 



(18) 



where we have assumed that the function $ tends to zero 
both at the center and the infinity. The imaginary part of 
equation (jlSl) gives 



= 9(m) 



7q\^\^jr*dr. 



(19) 



Since q is positive deflnite for < r < J? as seen in equation 
(p^, we have 9(/i) = 0, except the trivial case "I> = 0. That 
is, the eigenvalue /i must be real number. In a similar way, 
we introduce <I> — (ro — ^)/ to have another integral relation. 



0= f (^\f'\' +v\ff^jrUr. 



(20) 



The function within the integral is positive at least for < 
r < R. We therefore have the contradiction. 

If the eigenvalue is located within the domain ( p^ ) , the 
eigen-function has the singular point, say, r«, inside the 
star. The function is approximated by the delta-function 
as / ~ S{r — r»). The quasi-toroidal fluid velocity has the 
form U ~ 5{r — rt) from the equation (^. The function 
represents a steep resonance between the perturbation and 
the mean flow. The resonance may be more clear, if we con- 
sider so-called the Cowling approximation. In the Newtonian 
pulsation theory, gravitational perturbations are sometimes 
neglected in the oscillations. This trick gives good results 
for the spheroidal mode as well as the r-modes. The rela- 
tivistic Cowling approximation is given by ST^^.^ — with 
Sgap = 0. One component of the equations is reduced in the 
slow rotation case to 



{■CO - ii)U = 0. 



(21) 



It is clear that the solution of this equation is U ~ S{r — r^,), 
and that the range of eigenvalue is given by equation ( |l7| ) . 

Finally, we consider the Newtonian limit, in which zu 
Q. The frequency therefore corresponds to a single value as 



1(1 + 1) 



mCl. 



(22) 



This is the frequency of the r-mode oscillation measured by 
inertial frame. 



4 DISCUSSION 

In this paper, the r-mode oscillation is examined as the 
consistent first-order solution to the quasi-toroidal motion. 
The frequency forms a continuous spectrum. The oscilla- 
tion is caused by a certain resonance between the perturba- 
tion and the background rotating flow. The resonance con- 
dition is that the co-rotating frequency, [a — mn)e~'^^^ of 
the wave should be —2m/{l{l + 1)) times the angular veloc- 
ity, ■roe""''^ measured by ZAMO (zero-angular-momentum- 
observer). The angular velocity depends on the position of 
the local inertia frame due to the dragging effect. In this way, 
the r-mode oscillations in relativistic stars are much analo- 
gous to those in the differential rotating discs, although the 
angular velocity, Q is uniform. The mechanism works every- 
where within the rotating star, but the resultant frequency, 
a measured at infinity is not identical. This is the physical 
meaning of the continuous spectrum of the r-mode. 

The eigen-function of the Newtonian r-mode is not de- 
termined to first-order of the rotation, since any functions 
for the same fi satisfy the equation governing the oscillation. 



(^i- 0)17 = 0. 



(23) 



The modes are degenerated in this sense. In order to deter- 
mine the radial structure of the r-modes, calculation of the 
next order is necessary. The higher order corrections to the 
frequency will remove the degeneracy. As for the relativis- 
tic r-mode, the frequencies are distinguished corresponding 
to the resonance points. All the positions are on an equal 
footing to the first-order of the rotation. Therefore, the nor- 
mal frequency forms a continuous spectrum. We expect that 
some favored resonance points are selected as a result of the 
higher order corrections. That is, the axial part of the first 
order drives the density and pressure perturbations at the 
second order. The gravitational radiation may be also as- 
sociated with the density perturbations. The polar pertur- 
bations react on the frequency at the third order. The in- 
ternal structure will strongly affect the modes through the 
coupling. The relevant second-order rotational corrections 
~ 0{£^) like rotational deformation are of course necessary 
to solve the problem. The study beyond the flrst-order of 
the rotation is very important not only for the radial struc- 
ture, but also for the stability, although the calculations are 
significantly complicated. 

The frequency at the first order is a real number, and 
the mode represents standing ripple in the rotating flow. The 
wave will decay or grow due to the dissipation. The grav- 
itational radiation reaction and/or the viscosity cause the 
instability of spheroidal modes in the rotating star. Simi- 
lar instability may set in for the r-mode, according to the 
general argument (Friedman & Schutz 1978; Friedman & 
Morsink 1997). Recent numerical calculation suggested the 
instability of the r-mode (Andersson 1997). However, these 
works are not in agreement as for the growth rate, which is 
higher order consequence of e"{n > 2). 
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In conclusion, the second-order effect to the r-mode os- 
cillation in the relativistic star is complicated, but quite in- 
teresting problem. 
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